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We study the classical radiation emitted by free-falling charges in de Sitter spacetime 
coupled to different kinds of fields. Specifically we consider the cases of the electromagnetic 
field, linearized gravity and scalar fields with arbitrary mass and curvature coupling. Given 
an arbitrary set of such charges, there is a generic result for sufficiently late times which 
corresponds to each charge being surrounded by a field zone with negligible influence from 
^ ' the other charges. Furthermore, we explicitly find a static solution in the static patch 

' adapted to a charge (implying no energy loss by the charge) which can be regularly extended 

\^ i beyond the horizon to the full de Sitter spacetime, and show that any other solution decays 

at late times to this one. On the other hand, for non-conformal scalar fields the inertial 
O \ observers naturally associated with spatially flat coordinates will see a non-vanishing flux 

far from the horizon, which will fall off more slowly than the inverse square of the distance 
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^ij^' for sufficiently light fields (m^ + ^i? < 5i/^/4) and give rise to a total integrated flux that 

grows unboundedly with the radius. This can be qualitatively interpreted as a consequence 
, of a classical parametric amplification of the field generated by the charge due to the time- 

I dependent background spacetime. Most of these results do not hold for massless minimally 

' coupled scalar fields, whose special behavior is analyzed separately. 

(N 
cn 

\Q '. I. INTRODUCTION. 

O 

The study of interacting quantum field theories in de Sitter spacetime [iH^ has recently been 
the subject of an increasingly renewed interest [7-20]. There is controversy in the literature on 
^ , the significance of large quantum infrared effects and whether they can lead to an instability of de 

Sitter-invariant states and even a secular screening of the cosmological constant. In this respect, 
it is good to have solid answers to basic questions concerning field theory in de Sitter space. 

In this note we want to analyze in detail the possible classical radiation of scalar, electromagnetic 
and linearized gravitational fields by free-falling charges in de Sitter spacetime. The fact that a 
charge follows a geodesic does not necessarily imply the absence of radiation. This is clearly 
illustrated by the example of an electrically charged (or neutral) small mass orbiting around a 
central massive object in an asymptotically flat spacetime, which will emit electromagnetic (unless 
the orbiting mass is neutral) and gravitational radiation. In this case there is a natural set of 
asymptotic inertial observers with respect to which the orbiting mass describes a nontrivial motion, 
and the existence of radiation can be understood fairly naturally; the situation is, however, less 
clear for de Sitter. Note that the emission of radiation in the previous example will lead to a 
back-reaction force which will deviate the actual trajectory from the original geodesic motion. 
In contrast, the high degree of symmetry of de Sitter guarantees that the trajectory of a single 
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free- falling charge will not be altered even if it radiates. 

In flat space radiation is defined as a far-zone field which can be separated from the near-zone 
field associated with the sources [2]|, and falls off slowly enough so that it can propagate the 



total energy flux (integrated over the whole solid angle) over an arbitrarily long distance (in the 
absence of absorption or damping). Thus, in formulating the concept of radiation as emission of 
energy by the charges and subsequent energy propagation independently of them, local and global 
energy conservation arguments play a crucial role. These can be understood as a consequence of 
Minkowski spacetime possessing globally timelike Killing vectors. 

When considering radiation in curved spacetimes, one needs to confront two new conceptual 
difficulties. The first one is the absence of energy conservation. Even though the total stress tensor 
is covariantly conserved, one can no longer relate the change of the energy within a given region to 
the flux through the surface enclosing that region in a simple and direct way. This is still possible 
when there is a timelike Killing vector available (see Sec. Ill C"T]) . but this is not the case in general. 
It is, therefore, possible to introduce this kind of construction within a static patch of de Sitter, 
leading to a conserved Killing energy, but not globally. (In asymptotically flat spacetimes it is also 
possible to define consistently a notion of total energy and asymptotic energy flux in terms of the 
ADM and Bondi masses, and related concepts 22l].) 

The second new difficulty is the absence in general of a preferred frame of observers with respect 
to which one could naturally define the energy density and local energy fluxes (given any frame 
characterized by a tetrad {ca}, they are obtained by contracting the stress tensor with the tetrad 
appropriately). The importance and possible sensitivity to the choice of frame of observers can be 
illustrated with a simple example in flat space which shows that the flux can be zero for one set 
of observers, but nonzero for a different set. Consider an eternally uniformly accelerated electric 
charge in Minkowski space (with proper linear acceleration a). An inertial observer (placed at a 
distance much greater than 1/a from the charge) would see a nonzero flux. On the other hand, 
Rindlerian observers "comoving" with the charge would flnd a vanishing flux 23|. (The easiest 
way to see this is by realizing that in the Rindler frame one has a static charge which produces a 
static electric held 



23|. 



We will bear these issues in mind when studying the flelds generated by free-falling charges in de 
Sitter. As a consequence, we will separately analyze the energy loss by the source and the fall-off of 
the integrated flux at large radii (larger than the horizon), two aspects which are independent due 
to lack of global energy conservation. Furthermore, we will examine each case in several different 
frames. On the one hand, we will consider the inertial observers naturally associated with spatially 
flat and global coordinates, which are both particularly suitable for studying fluxes at radii much 
larger than the horizon. (Working with global coordinates can unmask global restrictions, such 
as the requirement for a vanishing total electric charge, which are absent otherwise.) On the 
other hand, we will also study the problem in terms of static coordinates, whose associated static 
observers (except for the central one) are non-inertial. Within the static patch covered by these 
coordinates one can deflne a conserved Killing energy and they are particularly useful to analyze 
the energy loss by the charge (which must have its origin in the current-field interaction term). 

Various aspects of classical fields and radiation of massless scalars in de Sitter space have been 



3 



studied in a number of papers 24427|| (certain quantum mechanical features of radiation emission 



in de Sitter have also been analyzed in Refs. [28|, l29|]). Here we will address the question of 
radiation of free-falling charges in quite some generality. We consider a broad class of fields (scalar 
with arbitrary mass and curvature coupling, electromagnetic and linearized gravity) and obtain 
the generic late-time field for an arbitrary set of free-falling charges. Furthermore, we examine 
the situation in several different frames to get a more complete picture, and pay especial attention 
to those aspects connected to radiation, namely, energy loss by the charge and fluxes at large 
distances. 

The paper is organized as follows. We start by considering the electromagnetic and conformal 
scalar fields in Sees. ITIl and lllli In both cases the calculation can be reduced to a calculation in flat 
space by using a suitable conformal transformation, and complete absence of energy flux is found. 
Next, the massless minimally coupled scalar field is considered in Sec. IIVI The charge is found to 
lose energy at a constant rate, but the total integrated flux falls off to zero at large radii. The 
case of general mass and curvature coupling is studied in Sec. |V] and solutions implying no energy 
loss by the charge are found for = + > 0. Nevertheless, for sufficiently light fields one 
has non-vanishing total flux at large radii. In Sec. IVII we analyze the case of linearized gravity. 
The generic late-time behavior for an arbitrary set of geodesies is described in Sec. IVIII Finally, 
we discuss our results in Sec. IVIIII Further details about coordinate systems in de Sitter and the 
regular extension of the static solution to the full spacetime are provided in a couple of appendices. 

We use the (+, +, +) sign convention of Ref. 30(] and work throughout the paper with natural 
units {h = c = 1) such that the Hubble constant equals one. 



II. ELECTROMAGNETIC FIELD 



We start with the calculation of the electromagnetic field produced by a free-falling charge (i.e. 
moving along a geodesic). Since any two geodesies in de Sitter can be related by an isometry 
transformation, we can consider without loss of generality the geodesic corresponding to 63 = = 
r = in the different coordinate systems (global, spatially flat and static) described in Appendix A. 
When working with global coordinates, we will need to include explicitly a second particle with 
opposite charge, which will be considered to move along the antipodal geodesic. Taking into 
account that Maxwell's theory is conformally invariant and de Sitter space is conformal to the 
Einstein static universe or to flat space (for the region covered by a spatially flat foliation), we 
expect no electromagnetic radiation to be emitted by a free-falling charge in de Sitter. We will see 
that explicitly in the next subsections. 

The action for the electromagnetic field coupled to a current is 

S[A^, J^'] = I d'x^ (^If^.F'^'' + A^J>^^ , (1) 

where F^^ = 2V[^74j^] = 2d[nA^j and the indices are raised with the inverse metric as usual. It 
leads to the following equation of motion: 

V^F^" = -^d^i^Fn = J": (2) 
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For a point-like charge with worldhne z^^{t) the current is given by 

J^(x) =q [dT ,5(4) - z"(t)) = ^(3) _ ^a^^^^ _ 

with ' = d/dr and where the first equahty is vahd for any parameterization of the worldhne, but the 
second one is only valid when r is the proper time and (locally) synchronous coordinates adapted 
to the worldhne are employed, so that gj^J is the spatial metric of the subspace orthogonal to the 
worldhne. The stress tensor for the electromagnetic field, which can be obtained by functionally 
differentiating the action for the free field, is 

T^u — FfiaFu ~^9pLvFal3F ^ ■ (4) 

As mentioned above, one can see that the action in eq. ([1]) is invariant under a conformal transfor- 
mation of the metric, g^j^u — )• g^j^u, while keeping (and hence F^y) fixed. Note, however, that 
the physical values of the electric and magnetic fields, which corresponds to contracting F^y with an 
orthonormal basis {ea} will transform as F^^Cae'^ — ?■ r2~^F^,yeae^ due to the rescaling required to 
keep the basis orthonormal. Similarly, although T^y would rescale with a factor ri~^, the physical 
values of the energy density, pressure and energy fiux will transform as T^pC^c^ — V ^T^^c^c^. 
Finally, from eq. ^ one can see that the current and its physical values transform respectively as 
follows: Jf" n-^Jf and J^ef, ^ ^.-■'Ji'ef,. 

A. Spatially flat coordinates 

As seen in eq. when parametrizing half of de Sitter space using spatially flat coordinates, 
it is conformal to half of Minkowski spacetime with conformal factor = l/rf. Therefore, the 
electromagnetic fleld for a free-falling charge with x* = const, [and hence dx^/dr] = (1,0,0,0)] will 
be given in the Coulomb gauge by^ 

^0=4^, A = and Fo.=q^, (5) 

where |x| = {dijX^x^Y^"^ and we only wrote the non-vanishing components of the field strength. The 
physical electric field, F^j = q{rf' / A'k\x\'^){x'^ /\x\) , points in the radial direction and its magnitude 
is inversely proportional to the square of the physical distance on spatial sections of constant time. 
Since there is no magnetic field in this frame, the energy fiux measured by the inertial observers 
comoving with the charge is zero, i.e. Tqi = 0. 



B. Global coordinates 

The situation is very similar when considering the global coordinates described in Appendix A.l. 
In this case, one can take advantage of the fact that the spacetime is manifestly conformal, with 



^ Note that the Lorentz gauge condition V ^lA^ = 0, satisfied by our solution in Minkowski spacetime, is not preserved 
in general under conformal transformations. 
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a conformal factor a^(A) = cos ^(A), to a region of the Einstein static universe with line element 
ds"^ = -d\^ + dVLl and -7r/2 < A < tt/2. 

Working in the Coulomb gauge, VjA* = = V^j4^, one can easily find the solution for a 
free-falling charge at ^3 = in the Einstein static universe: 

^0 = ;^cot(^3), Ai = Q and ^03= (6) 
47r 47rsm (t^3) 

where F^" can also be obtained directly from eq. ([2]). Because of conformal symmetry, the so- 
lutions have exactly the same form in global de Sitter coordinates (although no longer sat- 
isfies the Lorentz gauge condition in de Sitter space). The physical electric field is given by 
= g/47ra^(A) sin^(03), which is inversely proportional to the physical area of the 2-sphere cen- 
tered at the charge and contained on the 3-spheres corresponding to the constant-time spatial 
sections. Note also that the solution in eq. ([6|) implies the existence of a charge with opposite sign 
at the antipodal geodesic (^3 = vr). This can be understood as a consequence of Gauss's law [the 
component of eq. ([2|)], which implies that the total charge (i.e. the integral of the charge density 
over the compact spatial sections) should vanish. 

As in the case of spatially flat coordinates, the energy flux seen by the comoving inertial observers 
in the global frame vanishes (Toj = 0) due to the absence of a magnetic field. 



C. Static coordinates 



Working with the static coordinates defined in Appendix A. 3, one can also solve Maxwell's 
equations easily to obtain the following result in the Coulomb gauge (Vj^* = = V ^A^): 

^* = £:' ^^ = F*r = ^, (7) 

where the component Ft^ coincides in this case with physical value of the electric field corresponding 
to the component F^^. in an orthonormal basis. Once again the absence of a magnetic field implies 
a vanishing energy flux seen by the static observers (T^^ = 0). 

Note that the results for F^jj that we have obtained in the three coordinate systems are equiva- 
lent. This can be immediately seen by using eqs. ([70|) - (f75]) to show that the functional dependence 
of the only non-vanishing component of Ff^^ is the same in all cases, and taking into account that 
the transformation which relates the orthonormal bases associated with the different coordinate 
systems leaves the relevant element of the exterior basis invariant (the transformation in tangent 
and cotangent space is a boost along the direction of the electric field, which leaves it invariant). 
The solution in static coordinates also agrees with the result for the electromagnetic field of the de 
Sitter-Reissner-Nordstrom solution in the limit of very small mass and charge. 



1. Killing- energy conservation 



The existence of a timelike Killing vector in the static patch (whose tangent curves are 
parametrized by t) allows the definition of a conserved Killing energy constructed from the current 
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jfi = T^uC^, which is covariantly conserved (i.e. V^j/i =0). Applying Gauss's theorem to a cyUn- 
drical spacetime region with top and bottom boundaries corresponding to constant t hypersurfaces 
plus a third boundary generated by time translation, and differentiating with respect to time, one 
gets 

[ d'x^T,,en^^ = N [ d^x^T^A'k^ (8) 

where is a the unit vector normal to the hypersurface St, k^^ is the outward unit vector normal 
to its boundary 9St, and is the normalization factor of the Killing vector, defined by = Nn^. 
This equation displays the conservation of the Killing energy by relating the rate of change of 
the energy contained within the surface dTit and the total flux crossing that surface. Note that 
in principle one needs to consider the total stress tensor because it was necessary to consider a 
covariantly conserved stress tensor ( V^T^j^ = 0) when deriving eq. ([8]) . 

In particular one can apply this result to study the energy loss (emission) by the charge. In 
this paper we will focus on the contributions to the stress tensor from the field and current-field 
interaction term. (In the examples that we will be considering the trajectory of the free-falling 
charge will not be affected and the contribution of the free charge terms will not play a significant 
role.) As we will see, whenever one has a static solution for the field, the flux [the right-hand side 
of eq. dS])] vanishes, and so does the contribution of the free field to the left-hand side. From this 
we can conclude that the charge did not experience any energy loss. And indeed the only changes 
to the charge energy (if its trajectory is not altered) can come from the current-field interaction 
term, which will remain constant for a static solution. On the other hand, in Sec. IIVBI we will 
have a non-stationary solution with a non-vanishing flux for a massless minimally coupled scalar 
field, and this kind of energy conservation argument will provide an alternative way of obtaining 
the rate of change of the current-field interaction term. 



III. MASSLESS CONFORMALLY COUPLED SCALAR FIELD 



As we will see, the situation for a massless conformally coupled scalar field is very similar to the 
electromagnetic case. The action for a real scalar field with arbitrary mass m, curvature-coupling 
parameter ^ and coupled to an external source J{x) is given by 



S[4>, J,g^a 



(9) 



from which one obtains the following equation of motion, known as the Klein-Gordon equation: 

1 



(□ - - ^R) . 



-g 



(10) 



where R is the Ricci scalar. For a point-like charge with worldline z^{t) the current is given by 

^(4) _ ^«(^)) = ^(3) 



J{x) = q dr 



X — z 



j(3) 



[X — z 



(11) 
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where the first equahty is valid for any parameterization of the worldline whereas the second one 
is only valid when (locally) synchronous coordinates adapted to the worldline are employed. The 
stress tensor is obtained by functionally differentiating the action with respect to the metric 

Tf,^ = ^^(t>^u<P - \9nv (V"0V„</) + m202) + ^ ^g^^u - V^V^ + G^^) (12) 

For vanishing mass and = 1/6 (commonly known as the conformal coupling case), the action 
S[(j), J,g^u\ is invariant (up to a total divergence) with respect to conformal transformations of the 
metric, g^y — )• g^y = ^'^{x)g^y, together with the following rescalings of the field and the current: 
(j) ^ (j) = and J ^ J = Q~^J [note that the current for a point-like charge, given by eq. (jlip . 

does rescale appropriately] . It is then easy to see that the stress tensor transforms in a simple way 
under such conformal transformations: 

where we used the invariance of the action with respect to conformal transformations in the sec- 
ond equality. This simple transformation of the stress tensor under conformal rescalings will be 
exploited below to obtain the energy flux in de Sitter space from the result in Minkowski or the 
Einstein static universe. 



A. Spatially flat coordinates 

As mentioned above, for a massless conformally coupled field (m = 0, ^ = 1/6) the action is 
invariant under conformal transformations. Therefore, by taking = —l/rj one can immediately 
relate the field generated by a free-falling scalar charge in de Sitter when using spatially flat 
coordinates to the analogous situations in flat space, where the solution is given hy (j) = — g/47r|x|. 
The corresponding solution in de Sitter is then obtained by rescaling the field: 

^ «(-^) (14) 



47r \x 



Since the Minkowski solution is time independent, there is no energy flux associated with it (Tqi = 
0). Eq. ()13p then implies the absence of energy flux as seen by the comoving observers (x* = const.) 
in de Sitter spacetime with spatially flat coordinates. 

B. Global coordinates 

Similarly, the solution for a free-falling scalar charge at ^3 = in global coordinates can be 
straightforwardly obtained by solving the Klein-Gordon equation in the conformally related Ein- 
stein static universe, where the relevant equation to solve (if one looks for a static solution) is 



sm\e^)de?, V 003 J ^ 47rsin- 
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which follows from applying eq. (jlOp with m = and ^ = 1/6) to the Einstein static universe 
and taking into account that R = 6 when the radius of curvature is one. Introducing z = cos(03), 
eq. (fTSi) becomes 

d'^ch 86 q 5(z — 1) , . 

^ ' dz^ dz ^ A7tV2 ^ ' 

The homogeneous part of this equation corresponds to the hypergeometric differential equation 
with parameters a = h = \ and c = 3/2. Thus, the solution to eqs. (llSp and ()16p is 

. -ji, ,3/2; i±i) ^ -ji.-^ ^ (l - I) , (IT) 

which is the only solution of the homogeneous equation which is regular everywhere except for 
z = 1, where the source is located, and has the right behavior as z — )• 1 compatible with the right- 
hand side of eqs. (fT5]) and (fT6|) [namely, 4> ~^ —qj^T^Qj, as ^3 — )■ 0]. There is a second independent 
solution of the homogeneous equation, 



'47r2 



F(l,l;3/2; +Ffl,l;3/2;- 



47r sin(6'3) 



(18) 



which is not regular at z = — 1 either and actually corresponds to having also an identical source 
at the antipodal worldline (^3 = vr). This solution will play an important role when comparing 
with the solutions in the other subsections and in Sec. |Vl Other solutions of the homogeneous 
equation can be obtained by taking a linear combination of eqs. (jl7p and (jlSp and will correspond 
to a certain combination of sources at ^3 = and ^3 = tt. 

Note that in contrast with the electromagnetic case, where Gauss's law forced the total charge 
on the compact spatial section to be zero, even static solutions of the Klein-Gordon equation can 
have a non-vanishing total charge in the Einstein static universe due to the term proportional 
to the Ricci scalar and the field in eq. (fTUj) . as illustrated by the solution in eq. ([UD. In fact, 
by considering a continuous superposition of rotations on of that solution, one can generate a 
solution which corresponds to an arbitrary charge distribution on 5^. 

The solutions in de Sitter space can be immediately obtained from those in the Einstein static 
universe dividing by a(r) = cosh(r). For instance, the solution for two identical antipodal charges 
becomes 

^{T,Q^) = — Ji (19) 

47r cosh(r) sm(t'3) 

The static solutions in the Einstein static universe have obviously Toj = 0. Taking into account 
eq. ()13p . we conclude that there is no energy flux as seen by the comoving observers (with 0* = 
const.) in de Sitter global coordinates for the solution in eq. (fT9|) or any solution obtained by 
rescaling the static solutions mentioned above for the Einstein static universe. 
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C. Static coordinates 

Finally, one can work with the static coordinates of Appendix A. 3. If we consider a free-falling 
charge at r = 0, the static solutions will satisfy the following equation 



Unfortunately in this case it is not straightforward to find explicit solutions in terms of known 
special functions, but one can still study their existence and their main properties. This will be 
done in Sec. |V] for arbitrary values of ^ and m. As we will see, eq. (|20p has a one-parameter family 
of solutions, but only one of them is regular on the horizon (at r = 1). From the construction in 
Sec. Eland Appendix B, it is clear that such a solution should correspond to eq. p^ . which is 
regular everywhere except for the charge worldlines at ^3 = and 63 = tt. Using eq. ([73]) it can be 
written in static coordinates as 

^(*'-) = -£' (21) 

which can be easily checked to satisfy eq. (f20]) . Using eq. (fTOl) one can also see that it is equivalent 
to the solution in spatially flat coordinates given by eq. (jl4p . Being time- independent, this solution 
has zero flux as seen by the static observers with r = const, (i.e. Ttr = 0). 

IV. MASSLESS MINIMALLY COUPLED SCALAR FIELD 

The main features of the massless minimally coupled case are the absence of a stationary solution 
in the static patch which can be regularly extended to the full spacetime, and the existence, instead, 
of solutions which exhibit a constant non-vanishing energy flux as seen by static observers. Some 



of the results presented here have previously been obtained in Ref. 25(, where this problem was 
considered. 

A. Spatially flat coordinates 

In this case conformal symmetry is no longer available and we need to solve explicitly eq. (jlOp 
with m = ^ = 0, which becomes (after multiplying by \/—g) 

d^'^d,cP=^S^^Hx), (22) 

where we particularized the current in eq. (jlip to a free- falling point-like charge located at = 0. 
The retarded Green function for a massless minimally coupled scalar field in de Sitter spacetime 
is a well known result. An explicit derivation starting from eq. (|22p can be found for instance in 
Ref. [25,] and the result is 



11 1 

— S [r] — T]' — \x — x'\) + —9 [r] — rj' — \x — 



4tt\x — x'\ An 



e(r?-r?'), (23) 
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where one should remember that r/ is negative and future infinity corresponds to r/ = 0. Hence, we 
will have the following solution of eq. (j22p : 



d x' Gr,(7/, x; rj' , x' 



rj' Air 



+ log — 7— 



+ const. (24) 



where r/o is an arbitrary dimensionful parameter and the remaining constant term is actually 
divergent and reflects the logarithmic divergence exhibited by the integral with respect to rj' . Nev- 
ertheless, one can simply discard this constant term and it will still be a solution of eq. (j22p . Note 
that in contrast with Ref. 25|], our solution is valid for all times, although it becomes arbitrarily 
large at early times. This point and the possibility of "turning on" the charge at a finite initial 
time r]i, as done in Ref. [25!], will be further discussed in Sec. IIVBI 

Calculating the time and spatial derivatives of the solution found in eq. ()24p . 

doc^ = f f ^ + , (25) 



47r \\x\ rj — \x\ 

d(j) = ^ ( I \ (26) 

* 47r \ I^P \x\ {\x\ — rf) J ^ 

we can find the energy flux, whose components in an orthonormal basis are given by Tgj = 
ri'^{do(p){di(p). This flux is invariant under the isometry associated with the time translation in 
eq. (j65p . The first contribution to 4>, do4> and 9j(/> dominates for points well within the horizon 
{\x\ <C |r/|) and coincides with the solution for the massless conformally coupled case. Its contri- 
bution to Tqj is not zero, but falls off fast enough so that its total integrated flux vanishes in the 
limit of infinite radius. (In the conformal case Tg^ vanishes because there is an exact cancellation 
from the ■^V^Vjy^^ term.) For points outside the horizon (|x| > \rj\) the second terms on the 
right-hand side of eqs. ()25p -()26p must also be taken into account. Such a term dominates then the 
contribution to di(j) for large |x| and is inversely proportional to \x\. On the other hand, the two 
terms on the right right-hand side of eq. ()25p become comparable and partially cancel out, so that 
do(j) falls off like for large This means that the total integrated fiux tends to zero for 

large radii. [These points can be seen more directly by considering first a large distance expansion 
on the right-hand side of eq. (j25p .] 



B. Static coordinates 

In this subsection we will consider the static coordinates of Appendix A. 3, but it will be conve- 
nient to introduce the tortoise radial coordinate r*, which arises from requiring dr^/dr = (1 — r^)~^ 
and is related to the standard radial coordinate by r = tanhr*. In terms of this new coordinate, 
the line element becomes 

ds" = J" + tanh^ (r, ) d^l . (27) 
cosh (r*) 

We will consider the same free-falling source as in Sec. IIV A| whose trajectory is given in terms of 
these coordinates by x^{t) = (i, 0,0,0). The Klein-Gordon equation (including the source) for a 
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spherically symmetric field will be 

cosh(r,) \-tanh\n)d^<P + drJtanh'^{n)drA)] = -^^r^)- (28) 

47r 

If we are interested in finding static solutions, the equation further reduces to 

cosh2(r,) 92 + 2 coth(r,) = ^^t^ Hn) = (29) 

47r sinh (r*) 47rrj 

It is easy to verify that (j) = — (g/47r) [coth(r=i,) — r*] is a solution of this equation. Thus, it may 
naively seem that even for a massless minimally coupled field a free-falling charge in de Sitter can 
produce a stationary field, which implies Toj = 0, in its associated static patch. There are, however, 
several important additional aspects that should be taken into account. Let us start by pointing 
out that eq. (j28p has the following obvious time dependent solution 

^=--;^[coth(r*)-r* + at + 6], (30) 
47r 

for any real constants a and b. Next, one can use eq. (j70|) as well as eq. ([7T]) divided by eq. ([70]) to 
obtain the following relations between static and spatially flat coordinates: 

— T] = cosh(r*) e~*, 

|x| = sinh(r*)e~*. (31) 

Taking these into account, one can easily see that eqs. ([50]) and are equivalent (after discarding 
the divergent constant piece in the latter) if one takes o = 1 and b = log 1 770 1- But what is the 
explanation for the additional freedom in choosing a and b that one appears to have when working 
with static coordinates? This can be understood as follows. If we introduce the advanced and 
retarded times u = t + and u = t — r^, we have —r^: + at = v {a — l)/2 — u {a + l)/2, and we 
can clearly see that it will diverge on the future event horizon (which corresponds to u — )• 00 while 
keeping u finite) unless a = 1. Therefore, requiring that the solution is regular on the future event 
horizon and can be extended beyond it, completely fixes the existing freedom (up to the trivial 
additivity constant) and selects that solution already found in Sec. lIV Al which is not stationary in 
the static patch. There is no stationary solution that can be regularly extended beyond the event 
horizon, but as we will see in the next section, this is a peculiarity of fields with = m^ + ^i? = 0. 

Let us take the solution in eq. (jSOp with a = 1 and study the Killing energy and Killing energy 
fluxes along the lines presented in Sec. Ill C 1[ Even though the solution is not stationary, it gives a 
time independent result for the Tu component of the field stress tensor an its contribution to the 
first integral in eq. ([30|) . Taking this into account and using the same kind of arguments that led 
to the derivation of eq. we can infer the following two conclusions. First, the total energy flux 
going through surfaces of different radii should be the same and, second, it should equal the energy 
loss by the charge (coming from the current-field interaction term). The first point can indeed 
be checked if one takes N = l/cosh(r*), ^Jqqy, = tanh^(r*), Ttr, = — (g/47r)^ coth^(r=i,), = 
J^*, = cosh(r*) d.J:'^ and substitutes them into the right-hand side of eq. ([8]) to get a constant 
integrated flux (independent of the radius of the dY^t surface) equal to /Att. Furthermore, one 
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can also check that the decrease in the total Killing energy, dE/dt = —q^ /At:, that follows from 
eq. dS]) for such an energy flux does agree with the decrease of Killing energy associated with the 
current-field interaction term, which is given by qd(j){t,0)/dt = —q'^/An. 

Several additional remarks are in order. First, if we had included the standard kinetic term 
for the trajectory of the point-like charge with a mass m, we would have found (even allowing for 
a general trajectory of the charge) that the current-field interaction term can be interpreted as a 
contribution to the "effective mass" of the charge, given by meff (r) = m+q (/)(t(r), O) . Note that the 
field diverges at the origin and one needs to introduce a suitable regularization and renormalization 
procedure, which involves taking an infinite bare mass m that cancels the corresponding divergence 
(all this is necessary because we consider the strictly point-like particle limit). The important point 
is that within this context this procedure can be carried out in a time-independent manner, so that 
changes with time of the field will directly reflect changes of the effective mass; see Ref. 25| for 
further details, but keep in mind that we use the opposite sign for the current-field interaction 
term here. The second remark is that the solution in eq. pop with a = 1 diverges on the past 
event horizon (this is easily seen in terms of u, v coordinates, where the past horizon corresponds 
to n — )• — oo while keeping v finite), which is equivalent to the fact that the solution found in 
Sec. IIV Al diverged for r] — )• — oo. This issue was not present in Ref. 25|], where the charge was 
turned on at some finite time to, and it can be understood as follows: if the charge loses energy 
at a constant rate and one demands the mass to be bounded, it cannot be radiating for an infinite 
period of time. Hence, it was assumed in Ref. 25| that the charge was turned on at some initial 
time to, when the renormalized effective mass had a finite value ttiq. At later times, the effective 
mass would be given by mefr(i) = rno — q'^/ATT{t — to). We can see that the effective mass will 
vanish at a certain time and become negative after that. In Ref. [25] it was essentially concluded 
that the charge would simply disappear after that. In our opinion, however, the question is not 
so clear. Think of a free-falling charge which was also electrically charged. It would not radiate 
electromagnetically, but it would lose energy due to the emission of the minimally coupled field. 
At some point its mass would vanish, but it cannot simply disappear because of electric charge 
conservation. A nontrivial process is likely to happen at this point, but the actual details may be 
rather subtle and complex. We will return briefly to this question in the discussion section. 



C. Global coordinates 

Here we reanalyze the situation using global coordinates. We will only consider a single charge 
(this aspect will be further discussed below). The source on the right-hand side of the Klein-Gordon 
equation will then be given by (g/47r) (5(^3)/ sin^(03) cosh^(T). There is a nice expression for the 
retarded propagator in terms of the invariant interval for a pair of points {3l|: 

Gr(z) = i- [5(1 -z) + e{i - z)] X e{T - r), (32) 

where Z{x,x') = X'^{x)X^ {x') rjAB is the de Sitter-invariant interval associated with the pair 
of points X and x'. In terms of the coordinates introduced in Appendix A.l and for the 
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particular configuration with = considered here, we have Z{x^ x') = sinh(r) sinh(T') — 
cosh(r) cosh(r') cos(03). Integrating with the source, we get the fohowing intermediate expres- 
sion for the solution 



6yY + sinh(T) sinh(r') — cosh(T) cosh(r') cos(0; 
+ e(l + sinh(r) sinh(T') - cosh(r) cosh(T') cos(6l. 



(33) 



The argument of the delta and theta functions is zero when 



r ^ -lifccosh(r) sin(6l3) 

sinh(r) - cosh(r) cos(e3) ' ^ ' 

where the upper sign corresponds to the retarded time (T' < T) and the lower sign to the advanced 
one (T' > T). Moreover, there will be no solution when A < ^3 — 7r/2 (where A is the conformal time 
introduced in Appendix A.l), which corresponds to spacetime points that lie outside the causal 
future of the entire charge worldline. In the end one gets the following final result: 

" I "47rcosh(r) sm{e3) 

+ const. I 6'(A - 6I3 + 7r/2), (35) 



-1 + cosh(r) sin(6'3) 
sinh(r) — cosh(r) cos(6'3) 



where the same remarks given below eq. (|24p concerning the existence of a logarithmic divergence 
in the integration and the possibility of adding an arbitrary constant also apply here. In fact, it 
is considerably less cumbersome to evaluate the integral in eq. (|33p by changing to spatially flat 
coordinates (this can be done because the support of the integrand is entirely contained within the 
region covered by the spatially flat coordinates), and the calculation becomes then very similar to 
that of Sec. lIV Al Alternatively one can directly check that, when expressed in terms of spatially flat 
coordinates, eq. ([35]) reduces to eq. (|24l) . In order to show this equivalence one needs to make use of a 
couple of relations between spatially flat and global coordiantes: —l/ij = sinh(r) + cosh(r) cos(03), 
which follows from comparing the expressions for X° - in Eqs. §1^ and 1^, and -\x\/r] = 
cosh(r) sin(^3), which follows from comparing the expressions for (X^X^ij)^/^ in Eqs. §1^ and 
(j63p . In addition, one needs ot use the fact that the right-hand side of Eq. (j34p with the upper 
choice of sign is identically equivalent to (sinh(T) -|- cosh(T) cos(03))/(l -|- cosh(T) sin(03)). 

We close this subsection with several remarks. First, the solution in eq. ()35p diverges at A = 
^3 ~ "^72, which corresponds to the divergence at r/ — )• —00 or the past event horizon already 
found in Sees. IIV AlUVB] respectively. As discussed in Sec. IIVB) this can be cured by switching 
on the charge at a finite initial time (it can be done smoothly to guarantee that the Klein-Gordon 
equation is satisfied in the whole spacetime). In contrast with the electromagnetic or gravitational 
case, this is mathematically consistent because there is no symmetry that requires the charge to 
be conserved. Second, note that it was possible to find a solution for the full de Sitter spacetime 
even though we considered a single charge. While the Gauss constraint forbids this possibility in 
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the electromagnetic case and something similar happens also in the gravitational case (as briefly 
explained in Sec. IVI Cp . there is no such restriction for a scalar field. If we looked for time- 
independent solutions, ViV*(/) = would imply a vanishing total charge in the massless minimally 
coupled case, but the existence of extra terms involving time derivatives in the Klein-Gordon 
equation means that there is no restriction on the total charge when one allows for a general time 
dependence of the solution. On the other hand, demanding certain properties from the solution 
may sometimes imply the need to consider additional charges. For instance, we will see in the next 
section that having an identical antipodal charge makes it possible to have a stationary solution in 
the corresponding two static patches which can be regularly extended to the full de Sitter spacetime 
for -|- S,R > 0. However, as we will see, that is not possible for m? + ^i? = 0, and there was 
consequently no particular motivation for considering more than one charge here. 

V. SCALAR FIELD WITH ARBITRARY MASS AND CURVATURE COUPLING 

We have seen that, similarly to the electromagnetic case, one can have a free-falling charge 
coupled to a conformal scalar field in de Sitter spacetime which does not radiate. In contrast, for 
a massless minimally coupled scalar field it is not possible to have a stationary solution in the 
static patch which is regular on the horizon (and can be extended beyond it). A charge coupled 
to such a field keeps losing energy indefinitely at a constant rate. We would like to investigate 
how generic this situation is: is the absence of radiation a peculiarity of the conformal case closely 
tied to conformal symmetry, or is on the contrary the situation for the massless minimally coupled 
case exceptional? (The case of vanishing mass and minimal coupling is known to exhibit several 
other peculiarities due to large infrared effects.) We will show in this section that the generic 
situation is the absence of energy loss by the charge. On the other hand, we will see in Sec. IV Bl 
that for sufficiently light fields the integrated energy fiux seen by inertial observers at distances 
much larger than the horizon can grow unboundedly with the radius. The existence of these two 
seemingly contradictory results is a consequence of the lack of global energy conservation in de 
Sitter spacetime; see Sec. IVIlJ for a more detailed discussion. 

More specifically, given a scalar field with arbitrary mass m and curvature coupling ^, whose 
action and equation of motion are listed in eqs. ([U|)- pip . its dynamics for a fixed background 
spacetime with constant scalar curvature depends on the "effective mass" = m? + (,R. We will 
see that given two identical antipodal charges at ^3 = and ^3 = tt, there is always a globally 
regular solution for > which is stationary in the corresponding static patches. This solution 
exhibits a vanishing energy fiux as seen by the static observers. (Note that the antipodal source 
is entirely absent from the spacetime region covered by the spatially fiat coordinates or a single 
static patch.) 

It is worth pointing out that having > does not necessarily imply that the fields cannot 
exhibit massless behavior: it is not trivial to define unambiguously the notion of mass when the 
Compton wavelength is comparable or larger than the radius of curvature. In fact, the case that 
shares some key properties with massless fields in flat space (namely, a retarded propagator with 
support entirely contained on the light-cone for conformally flat spacetimes such as de Sitter) is 
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the conformal case, which corresponds to = 2 (where we have used that i? = 12 in units with 
H = l). 



A. Static coordinates 

The Klein-Gordon equation in static coordinates with a point-hke source at r = is 

' '''^ii(rH^-r^m^^-M^,^m (36) 



where A2 is the Laplacian on the 2-sphere, covered with the coordinates {0i,^2}- If we look for 
static and spherically symmetric solutions centered at r = 0, the equation reduces to 

The homogeneous part is a second-order ordinary differential equation with regular singular points 
at r = 0, r = 1 and r = co. Thus, Puchs's theorem guarantees that one can use the Probenius 
method to construct two independent solutions as a couple of convergent series. The two solutions 
around r = can be chosen to be 

Mr) = ^a„r" = l + — r2 + ... , (38) 

n=0 



Mr 

r 

n=0 



- bur"" + c ln(r) (r) = - + f ) r + • • • . (39) 



The first one is regular at r = since {d(f)i/ dr){<d) = (otherwise it would have a spike at the origin 
of spherical coordinates and would not be differentiable there) and is, in fact, an even function of r. 
It is, therefore, a solution of the homogeneous equation with no source at the origin. On the other 
hand, the second solution (in this case c = 0) corresponds to a Dirac delta source at the origin. 
With the right normalization, namely (f){r) = —{q/4Tr)(f)2{r), it is a solution of eq. (j37p including 
the point-like source. 

We are interested in solutions which are regular on the horizon. So let us study the series 
expansion for the solutions around r = 1 and then discuss how they match the solutions found 
above. In this case the two independent solutions can be chosen as 



^a„(r-l)" = l-— (r-l) + ... , (40) 



2 

n=0 



Mr) = J^6„(r-l)" + ln(r-l),^i(r-l) 

n=l 

= ln(r - 1) - — (r - 1) ln(r - 1) + f — - 4 j (r - 1) + • • • . (41) 

The second solution is singular at r = 1. Hence, we should consider a solution proportional to 
01 (r), which is finite and regular at r = 1. In fact, it can be extended smoothly to the full de 
Sitter space, as discussed in detail in Appendix B. In general, this solution will match a linear 
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combination of the solutions found above: C</>i(r) = A(j)i{r) + B(j)2{r)- Taking B = —{q/ATi)^ so 
that it is a solution of eq. (I37p including the delta source, completely fixes A and C . Note that one 
should exclude the possibility that B = because this would imply that a solution of the equation 
with the source at r = would always pick a contribution proportional to 4'2{t) and would have 
singular behavior at r = 1. In Appendix B we prove that this possibility can always be excluded 
for > 0. 

From the previous paragraph we can conclude that for > one can always find a static 
solution (j)s{r) associated with a free- falling charge in de Sitter which is regular on the horizon and 
can be extended to the full spacetime, as discussed in more detail in Appendix B and the next two 
subsections.^ Since it is time-independent in static coordinates, this solution exhibits zero energy 
flux as seen by the static observers {Tu = 0). 

We close this subsection by briefly mentioning what things would have looked like if we had 
used the tortoise coordinate introduced in Sec. IIVBI The Klein-Gordon equation would be the 
same as eqs. (f28]l and (|29|) with an extra M'^tf) term on the left-hand side. Close to the source at 
r = 0, r and are almost the same and the solutions in terms of r* have the same form as in 
eqs. (p8|) - (j39|) . On the other hand, close to the horizon, which corresponds to — )• oo, the Klein- 
Gordon equation becomes approximately d'^(f)/dr1 = 0, whose solutions are a linear combination 
of (t>i{rtf) = 1 and 02(''*) = These correspond to the solutions in eqs. ([I0|) - ()1T]) . and from 

this point on the argument would proceed exactly as above. In addition, one could study how 
perturbations with respect to this solution which are initially regular on the horizon decay at late 
times. This amounts to imposing the boundary condition of vanishing flux from r^, — )• oo and is 
analogous to the quasinormal-mode analysis of perturbations outside a Schwarzschild black hole. 
We will, however, study the decay of perturbations using spatially flat and global coordinates. 

B. Spatially flat coordinates 

In spatially flat coordinates the Klein-Gordon equation for a field with an effective mass 
and a point-like source at x* = is 

-|:^-3i.|^..-.«^-MV^,.-'(.). (42) 

Given a solution (t)s{j') in static coordinates like the one found in Sec. IV Al (and its extension to 
r > 1), one can explicitly check that (/>s(e'^|x|) is a solution of eq. (|12]) . as it should be. Furthermore, 
its large-distance behavior can be obtained from eqs. ([79|) - (f80]) and is given by 

^,{t,x') ^ {e'\x\y^~ , (43) 

with A_ = 3/2 - ^9/4 -M2. 



^ In fact, one can consider r > 1 in eq. p7p to study the extension of the solution to the remaining two non-static 
patches and introduce an analogous expansion in terms of 1/r around r = oo, which is also a regular singular 
point. In that way, one obtains the same kind of late-time decay that will be described in the next two subsections, 
as well as the large-distance behavior for scales much larger than the horizon; see Appendix B.2 for further details. 
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Interestingly enough, one can establish the generic late-time behavior of an arbitrary solution 
and show that it tends at sufficiently late times to (/)s(e'^|x|). This follows from the fact that for 
> the homogeneous solutions of eq. (1421) decay exponentially at late times. Indeed, if we 
Fourier transform the comoving spatial coordinates, the homogenous part of the Klein-Gordon 
equation becomes 



After a particular mode leaves the horizon (i.e. when \k\e'^ > 1) one can neglect the contribution 
from the term proportional to k'^ in order to study the evolution of that mode. The solution at 
sufficiently late times is, thus, a linear combination of two exponentially decaying terms, exp(— A±r) 
with A± = 3/2 lb y^9/4 — M^, which always have a positive real part for > 0. In contrast, for 
= we have A_ = and one of the two solutions becomes constant outside the horizon. This 
can give rise to large infrared effects as modes keep leaving the horizon and is closely related to the 
peculiar behavior that we found for the massless minimally coupled case, including the absence of 
a regular static solution in the static patch. 

Taking into account the explanations above, one can see that the approximate result for (j)s{T, x*) 
in eq. ()43p can be employed to calculate the energy flux at sufficiently late times and sufficiently 
large physical distances from the source. Its time and spatial derivatives are given by 



d<f>s 



~ (-A_)(e-|f|)''", (45) 

The energy flux can then be obtained from eq. (jl2p . Let us consider ffist the minimally coupled 
case (.^ = 0). Only the first term on the right-hand side of eq. (|12p will contribute to the energy 
flux. At sufficiently late times and sufficiently far from the source, the result for the radial flux 
when working with an orthonormal basis is the following: 

T^,^{X.f[e^\x\) . (47) 

For positive but sufficiently small masses (0 < <^ 1), A_ will be positive but very small. 
This means that the amplitude of the stress tensor will be quadratically suppressed, but, more 
importantly, it will fall off almost as slowly as the inverse of the physical distance. Therefore, the 
total flux integrated over a sphere of radius R will grow almost proportionally to R for sufficiently 
large radii. This is a remarkable result whose interpretation and discussion is provided in Sec. IVIIIi 
In fact, in order to have a non- vanishing integrated flux at large radii it is sufficient to have 
< 5/4, which corresponds to A_ < 1/2 (M^ = 5/4 would give rise to a finite integrated flux 
in the limit of an infinite radius, whereas for smaller masses it would grow unboundedly with 
the radius). For non-minimal coupling (^ ^ 0) there is an additional contribution from the term 
(.^fi^u'P'^ in eq. ([12]), but the conclusions remain unchanged. 
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C. 



Global coordinates 



In global coordinates the Klein-Gordon equation has the following form: 



q 5{e^) q 6(93 - tt) 
47r 6*2 47r (6*3 - 7r)2 



(48) 



ar2 " cosh2 T 



2 ' 



where A3 is the Laplacian of the 3-sphere, parametrized by {^1,^2)^3}) and we included a free- 
falling source at ^3 = as well as an antipodal one at ^3 = vr, which is necessary in order to have 
a static solution in the static patch which is globally regular. This is discussed in Appendix B, 
where we also study its late-time and large-distance behavior, which is generically found to be 



and is valid for cosh(T) sin(03) ^ 1. Note that at sufficiently late times, as the sources get farther 
and farther way, one will have non-negligible values of the field in two disconnected regions around 
each one of the sources with very similar form to that found above for spatially flat coordinates. 
As will be discussed in Sec. IVIIl this can actually be generalized to an arbitrary set of geodesies. 

Similarly to the case of spatially flat coordinates one can also establish that the generic late- 
time behavior of an arbitrary solution of eq. ()48p tends at sufficiently late times to (/>s(r, 63). The 
homogeneous solutions of eq. (j48p can be decomposed in terms of spherical harmonics yLAii^) for 
the 3-sphere with time-dependent coefficients which are proportional to (coshT)^^/^ times a linear 
combination of the Legendre functions P^^-|^^2'^^(tanhT) and Q^'^y^'^ {tanhT), and decay at late 
times in the same way already found above for spatially flat coordinates. 

The worldline of the antipodal source lies entirely outside the spacetime region covered by the 
spatially flat coordinates considered in Sec. IV B[ We close this subsection briefly discussing its 
implications for the results found there. When solving eq. ([I2]) . the effects of the antipodal source 
would enter through the initial value data specified on one of the flat spatial sections. At sufficiently 
late times any solution of eq. ()42p will decay to (j)s[e^\x\), which has no information left on the 
antipodal source. There is actually a simple interpretation of this point. Suppose that we start 
with 03 < TT for the second source and then take the limit 63 ^ it. Doing so while considering a 
flat spatial section with constant r, means that the location of the second source will be taken to 
an infinite physical distance on that section. Since the field generated by the second source will 
fall off at large distances in the same way as given by eq. (|43p . it will not contribute in the limit of 
infinite separation (at least for > 0). That explains why there is no contribution to 0s(e'^|x|) 
from the antipodal source. 



Gravitational radiation can be studied in the context of linearized gravity, regarded in this case 
as small metric perturbations around a de Sitter background geometry: g^i, = g^^^ + h^^. One 
can then evaluate the usual geometrical objects for the perturbed metric, expand in terms of the 




(49) 



VI. LINEARIZED GRAVITY 
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metric perturbation and retain the terms strictly linear in It is convenient to express them as 



R^^ = ^(VaV,/i^ + v„vx- v.v^/i- v"v,V)> (50) 

where h = g'^^hai3, is the covariant derivative associated with the background metric and 
we will be using the background metric to raise and lower indices unless explicitly stated otherwise. 
An additional conceptual difficulty in this case is the existence of a gauge symmetry corresponding 
to transformations under diffeomorphisms. One can employ the background field method, as we 
did above, which preserves general covariance in terms of the background geometry. Nevertheless, 
the gauge freedom under transformations of the form h^i, — )• h^i, + 2V(^.^,^) needs to be fixed (it 
can be interpreted as the ambiguity in splitting the background and the perturbation). Even if one 
does so by imposing a gauge-fixing condition compatible with the background general covariance 
(background field gauge), in general the results will still depend on the particular gauge fixing. 
Only truly gauge-invariant physical observables will be completely independent. 



In order to study gravitational radiation, we will follow Ref. 3^ and compute Gu^'^\ which 
results from evaluating the Einstein tensor for the perturbed metric ^^j^ and keeping the terms 
strictly quadratic in h^^, in a particular gauge. This object suffers from the gauge dependence 
mentioned above: while being generally covariant with respect to the background geometry, it 
will depend on the particular gauge which was chosen for h^i,. (In asymptotically fiat spacetimes 
it can still be employed to calculate gauge-invariant quantities such as the ADM mass, the total 
angular momentum, or the asymptotic energy fiux, but that is not the case here.) Notwithstanding, 
we will follow this common approach and leave deeper scrutiny for future work. In this respect 
it will likely be useful to consider the linearized Weyl tensor C^Jpa, which is gauge invariant 
because its unperturbed counterpart vanishes and which will completely and naturally characterize 
the perturbed geometry at linear order. ^ It might then be worth exploring the use of the Bel- 



Robinson super-energy tensor 33|, which involves terms quadratic in the Weyl, as a possible way 



of characterizing linearized gravitational radiation in a gauge invariant fashion. 

A. Spatially flat coordinates 

The linearized Einstein equation for metric perturbations around a de Sitter background can 
be written as 

G;:« = 8^Gr;w, (5i) 



^ Together with the hnearized Weyl tensor, the linearized Riemann tensor is entirely determined by the linearized 
Ricci tensor, which can be cast in a gauge invariant form by raising one of the indices (the background quantity is 
then proportional to 5^, whose Lie derivative with respect to any vector vanishes). When one imposes the Einstein 
equation, this quantity vanishes outside the sources: 
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where the gravitational coupHng constant G can be regarded as the expansion parameter and the 
hnearized Einstein tensor is given by 

g;:^^) = g'^R^'} - h-^^R^^i - Is:(r^J'^ - h-^r^^I), (52) 



with given by eq. ([5T]) and R^^u = "^Qiiv being the Ricci tensor of the de Sitter background. 
Note that when perturbing the Einstein equation with one of the indices raised, the contribution 
of the cosmological constant term is entirely contained in the background equation and it does not 
appear in eq. ([5T]) for the linear perturbations. As for the background stress tensor T^^'^^ = g^pT^^-^ 
of a free-falling point-like mass, one can take 

= / ^^7^ '^^'^ - '"^^^ = ^^'^ (^'^ - ' 

where we have used the same notation as in eq. ^ and = g^^yZ^z'^ . 

In order to find an explicit solution in spatially flat coordinates, we will make use of the results 
derived in Ref. 



27|. Defining ijjy = h'i — {h/2)5u and imposing the gauge condition V ^ijju 



— {2/r])ip^, a closed set of equations governing the dynamics of x^J.v = Vuafp'jl were obtained in 
that reference. One can study the solutions of the inhomogeneous and homogeneous equation. 
The homogenous equation satisfied by the two physical polarizations of the gravitational waves 
correspond to the Klein-Gordon equation for a massless minimally coupled scalar field. Here we 
will, however, focus on the solutions of the inhomogeneous equation. For a point-like mass located 
at X* = 0, the only non-vanishing component that results from evaluating eq. (j53p is 



(0) 

00 



m(— r/) (J*-^-* (x*), where we took into account that z^ = e '^6q = —tiSq in the {?/, x*} coordinates 
that we are employing. The nontrivial part of the Einstein equations is then 

•9^X00 + -^oXoo - ^Xoo = -leyrGTo^Q^ 

or, equivalently, d'^i — ] =-167rG^, (54) 
\ 1] J V 

where d'^ = rj'^'^d^di, is just the fiat space D'Alembertian. A particular solution of eq. ()54p is the 
following: 

Xoo 4Gm AGmr] 

= -FT ^ ^0 = -Xoo = p^, (55) 

T] \X\ \X\ 

which implies 

h'o = -^-^, K = ^^l (56) 



Following Ref. 32 1 we will consider = (l/SirG) Gu^'^^ as the stress tensor of the linearized 
gravitational field. It can be directly obtained by evaluating the exact Einstein tensor for the 
perturbed metric g^i, [using the result in eq. (|56]) ] and keeping only the terms quadratic in the 
gravitational coupling constant G. In this way we find 

which gives a total integrated flux over the sphere that falls off to zero at large radii. 
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B. Static coordinates 

The perturbative solution for a point-like particle at r = in static coordinates can be obtained 
by expanding perturbatively in G the exact Schwarzschild-de Sitter solution 

ds^ = -{l-'-^- r^) dt^ + + dnl (58) 

Keeping only terms linear in G one gets the following result for the linear metric perturbations: 

2Gm 2Gm 

Now we should use this metric perturbation to construct g^u and calculate = {l/SnG) Gu^'^^ ■ 
Alternatively one can consider the exact Einstein tensor for the full Schwarzschild-de Sitter solution 
and keep those terms quadratic in G. In fact, since the Schwarzschild-de Sitter solution satisfies 
the equation G'j^ = 3 5^ (remember that A = 3 in our units), there is no dependence on G left 
when evaluating the full Einstein tensor outside the sources. This will, however, include a possible 
contribution to the full Einstein tensor from terms quadratic in G in the exact Schwarzschild-de 
Sitter metric, which needs to be subtracted. Therefore, one needs to take /i^^J = 4G^m^/r^(l — 
r^)^ and use eq. ([50]) to calculate its contribution to the perturbed Ricci tensor, and from it 
get immediately the Einstein tensor contribution to be subtracted. Fortunately, one can easily see 
that the Oi components will vanish and we can conclude that = 0, as one intuitively expects 
from the existence of a static solution. 



C. Global coordinates 

Similarly to the electromagnetic case, one also needs to consider at least two sources in the case 
of linearized gravity. Analogously to the requirement of a vanishing total electric charge implied by 
the Gauss constraint, the constraints from the Einstein equations impose the following restriction 
on the lowest-order stress tensor of the sources: 

d'x^TlS)n^C'' = 0, (60) 

St 

given a compact Cauchy hypersurface Sj with normal unit vector n^, and a Killing vector S,'^ . (It is 
crucial to have a compact Cauchy hypersurface; that is why this issue arises in full de Sitter space, 
whose Cauchy hypersurfaces have an topology, but not in the subregion covered by spatially flat 



coordinates.) This restriction is closely related to the phenomenon of linearization instability [35 1 
in this context and is sometimes referred to as "linearization instability constraint" . In particular, 
for the timelike Killing vectors in de Sitter this condition is not satisfied by a single particle, but it is 
fulfilled by two antipodal free-falling particles with equal mass. In fact, the simplest way to obtain 
the linearized gravity solution for this case is to consider the small mass limit of the Schwarzschild-de 
Sitter solution, restrict oneself to the region outside the black-hole horizon (including the extension 
beyond the cosmological horizon) and treat the mass perturbatively there. In this way, we can 
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directly obtain the metric perturbation in terms of the static background coordinates from eq. (jSSp . 
The result (obtained in the previous subsection) will cover the two antipodal static patches of de 
Sitter, but it can be regularly extended to the full de Sitter spacetime by considering the extension 
of the static coordinates to r > 1 as described in Appendix A. 3, and even changing to global 
background coordinates using eqs. ([73|) - ([75|) . 

VII. LATE-TIME RESULTS FOR AN ARBITRARY SET OF GEODESICS 

At sufficiently late times a general set of geodesies in de Sitter spacetime will become arbitrarily 
far apart. We will argue that in most cases the field generated by a general set of free- falling charges 
moving along such geodesies will exhibit a generic form at sufficiently late times. 

A. Scalar field 

The basic idea is simple. We saw in Sec. IV Bl that for > the field generated by a single 
free-falling charge tends at late times to a generic solution 0s(t, x*), which in turn falls off at 
large distances. The generalization to multiple charges is straightforward: since we are considering 
linear fields, one can simply apply the superposition principle. (One would need external forces to 
compensate the mutual forces and guarantee that the charges continue to follow geodesies, at least 
until the mutual forces can be neglected, which is the late-time regime that we are really interested 
in here.) Furthermore, at sufficiently late times so that all the geodesies are sufficiently far apart 
[and taking into account the fall-off of (?i)s(T, x*) at large distances], one can focus entirely on the 
field generated by the closest charge and neglect the effect from all the other charges. 

Therefore, the generic late-time field configuration would correspond to effectively non- 
overlapping field zones surrounding each one of the charges, where one could employ adapted 
spatially flat and even static coordinates (for points within the horizon of the charge) such that 
the results of Sees. IV AIIV Bl would directly apply. This result, however, does not hold for = 0, 
as we found in Sec. llVl when studying explicitly the massless minimally coupled case. On the other 
hand, for < <C 1 the conclusion will be true, but one will need to wait for very long times 
to make sure that the solution associated to each one of the charges decays to (/>s(r, x*) and that 
the different charges are separate enough. This is because A_, which controls the decay rate of 
any homogenous solution of the Klein-Gordon equation as well as the fall-off of (/>s(t, x*) with the 
distance, is very small in that case. 

B. Electromagnetic field 

The situation for the electromagnetic field is similar to that for scalar fields with > 0. We 
saw in Sec. |II] that there is a static solution in the static patch which can be regularly extended 
to the full spacetime. The remaining question is then whether the solutions of the homogenous 
Maxwell equation decay at late times analogously to the solutions of the Klein-Gordon equation. 
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This can be easily investigated by taking advantage of the invariance of Maxwell's equation under 
conformal transformations. For instance, in spatially flat coordinates the solutions for F^i, are the 
same as in flat space, i.e. free electromagnetic waves. However, the physical electromagnetic field 
(referred to an orthonormal basis) should be rescaled by the inverse square of the scale factor: 
Pfiu = rfP^y. Physically this corresponds to the amplitude of the electromagnetic waves being 
redshifted away as they propagate in the expanding spacetime, and results in any initial solution 
decaying to that of Sec. HI] at sufficiently late times, or the corresponding generalization for an 
arbitrary set of free-falling charges analogous to that described for scalar fields in the previous 
subsection. 

C. Linearized gravity 

Things are slightly more subtle for linearized gravity. As seen in Sec. IVI C\ there is a static 
solution which can be regularly extended to the full spacetime. However, whether other solutions 
generically decay to it at late times is less obvious. After an appropriate gauge-fixing and field 
redefinition, one can see that the homogeneous equation satisfied by the two physical polarizations 
of the linear metric perturbations are that of a massless minimally coupled scalar field [3^. This 
would seem to imply that given a non-vanishing initial condition, it will not decay at late times. 
Nevertheless, one should consider gauge invariant quantities which naturally describe the local 
geometrical properties. The linearized Weyl tensor, C^Jpa, is one such object which fully charac- 
terizes the perturbed geometry to linear oder. One can check that the linearized Weyl tensor for 
the homogeneous solutions of the linearized Einstein equation does decay at late times, so that the 
static solution and its extension indeed correspond to an attractor geometry at late times. 

VIII. DISCUSSION 

We have analyzed the classical fields generated by free-falling charges in de Sitter spacetime 
coupled to different kinds of fields. We were able to provide explicit results for the electromagnetic 
case, linearized gravity, and scalar fields with a wide range of values for the mass and the curvature 
coupling (M^ = m? + > 0). Furthermore, we have argued that at sufficiently late times one 
should have a generic result for an arbitrary set of free-falling charges in all cases (other than a 
scalar field with = 0), which would correspond to each charge being surrounded by a field zone 
with negligible influence from the other charges and where the results obtained for a single charge 
would apply. 

Except for scalar fields with = we have generically found "absence of radiation" in the 
static patch associated with each free-falling charge. More precisely, we have obtained a static 
solution with vanishing energy flux with respect to this frame. This implies conservation of the 
Killing energy and, together with the static character of the fleld, can be exploited to argue that 
there was no energy loss (emission) by the source. As discussed in Sec. IV Al and Appendix B, 
having a static solution in the static patch is simple, the non-trivial aspect is whether it is regular 
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on the horizon so that it can be extended to the fuh de Sitter spacetime. 

Furthermore, we have also studied the energy flux at distances much larger than the horizon as 
seen by the natural set of inertial observers in spatially flat coordinates, and found a remarkable 
result. For sufficiently light fields (with < 5i?^/4) the energy fiux does not fall off fast enough 
(slower than the inverse square of the distance) and the total flux integrated over a sphere of 
radius R becomes arbitrarily large for large radii. This phenomenon can, in fact, be qualitatively 
understood as the result of a classical parametric amplification process. In contrast with the Killing 
energy in the static patch, there is no energy conservation in the natural frame associated with the 
spatially fiat coordinates (and in any case there is no frame where energy is conserved for scales 
larger than the horizon), and one can interpret that the field generated by the free- falling charge is 
being amplified by the time-dependent background spacetime, analogously to what would happen 
to a harmonic oscillator with a time dependent frequency (e.g. a pendulum with a time-dependent 
length) . 

Note also that the unbounded growth of the integrated fiux with the radius does not necessarily 
entail dramatic effects for two reasons. First, the stress tensor itself is small (the flux component 
is quadratically suppressed by the mass of the light field), so that its contribution to the back- 
reaction on the spacetime geometry will be locally small (although a more detailed analysis would 
be required to make sure that there are no significant secular effects). Second, while this kind of 
phenomena would face an insurmountable tension with energy conservation in flat spacetime, it 
becomes much easier to avoid contradictions in de Sitter, where energy is not globally conserved. 

On the other hand, for scalar fields with = there is no stationary solution in the static 
patch which can be regularly extended to the full spacetime. Instead we explicitly found for the 
massless minimally coupled case that the charge loses energy at a constant rate. As discussed in 
Sec. IIVBI this means that its effective mass will vanish and start to become negative within a finite 
time. In Ref. [25[ it was assumed that the charge would simply disappear at that point. However, as 



illustrated in Sec. lIVBj with the example of a particle which is also electrically charged, we believe 
that the question can be more subtle and should be studied within an appropriate microscopic 
model of the charges. It seems clear that a negative mass signals some kind of instability, but its 
precise nature is less obvious. For instance, if one tried to model the microscopic constituents of 
the charge by massless fermions with a Yukawa coupling to the field, it is known that the theory 
exhibits an instability unless there is a minimum amount of nonlinear self-coupling for the scalar 
fleld. It might be that for certain parameter ranges the theory itself is unstable in de Sitter space, 
or that the effects of the non-linear terms in the scalar field potential need to be taken into account 
when considering the late-time evolution of the effective mass. 



Acknowledgments 

E.T.A. would like to acknowledge discussions with O. Kancheli and Y. Khriplovich. E.T.A. 
would also like to thank the Albert-Einstein-Institut, where most of this work was carried out, 
for hospitality. The work of E.T.A. was done under the partial financial support of grants for the 
Leading Scientific Schools NSh-6260.2010.2 and RFBR 08-02-00661-a. 



25 



APPENDIX A: COORDINATE SYSTEMS IN DE SITTER SPACE 



Four-dimensional de Sitter space with unit radius of curvature {H = 1) can be defined as tlie 
liyperboloid given by t/abX^^^ = 1 and embedded in tlie five-dimensional Minkowski spacetime 
with metric tjab = diag(— 1, 1, 1, 1, 1). It is invariant under the ten-dimensional 0(4, 1) group of 
isometries, which corresponds to the Lorentz group of the embedding spacetime. 



1. Global coordinates 



One can introduce a coordinate system {T, 6i} which covers the whole hyperboloid and corre- 
sponds to the following embedding: 

X° = sinh(T), X'' = -cosh(r)a;'', a = 1,2, 3, 4, 
u!^ = sin(03) sin(6'2) sin(0i), 

= sin(03) sin(6'2) cos(0i), 
00^ = sm{6s) cos(02), 

= cos(03), (61) 
where Sah^j""^^ = 1 and defines a 3-sphere. The line element in these coordinates is 

ds^ = -clT^ + cosh2(T) (inl = 1— l-dX"^ + dnl] , (62) 

^ ^ cos2(A) L ^ ^ 

with dQl = dOl + sin^(^3) [dO^ + sin^(6'2) dOl] . The last equahty in eq. ([62]) shows exphcitly that 
de Sitter space is conformal to a region of the Einstein static universe characterized by the finite 
range of the conformal time, — 7r/2 < A < ti/2. Six of the ten independent Killing vectors are 
manifest in these global coordinates, namely, those associated with the 0(4) isometry group of the 
3-spheres that correspond to the spatial sections with T = const. 

The worldlines with 6i = const, are timelike geodesies with proper time T and can be regarded 
as inertial, or free-falling, observers. (Note incidentally that any pair of geodesies can be related 
by a de Sitter isometry transformation.) Throughout the text we often consider the particular 
geodesic defined by ^3 = 0, and the antipodal one, with 6*3 = vr. 



2. Spatially flat coordinates 



The following alternative parameterization covers half of the hyperboloid, corresponding to 
— X'^ > (the other half can be covered by an analogous parameterization which changes the 
sign of the X^ and X^ coordinates): 



X' 



2r] 2 



X' = - — , i = 1,2,3 
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with r] G (—00,0). The hne element can then be written as 

ds^ = -dr^ + e^'^dijdx'-dx^ = \{-drf + 6ijdx^dx^), (64) 

where the comoving time r is related to the conformal time r] hy rj = — exp(— r). In terms of 
the conformal coordinates {??, x*} it is obvious that this half of de Sitter is conformal to half of 
Minkowski space with ij G (— oo,0). Furthermore, seven of the isometry generators of de Sitter 
space are manifest in this coordinate system. In addition to the six generators of the Euclidean 
group E(3) for the spatial sections with r] = const., the metric in eq. (|63p is manifestly invariant 
under the following transformation: 



rj ^ e "'rj or, equivalently, r — )• r + a. (65) 

The worldlines with = const, are timelike geodesies with proper time r, and here the particu- 
lar geodesic corresponding to ^3 = in the previous subsection is given in conformal coordinates by 
x'^(r) = (— e""^, 0, 0, 0). In contrast, the antipodal geodesic located at ^3 = vr is entirely contained 
in the other half of de Sitter space not covered by this parameterization. 

3. Static coordinates 

While de Sitter spacetime has no global timelike Killing vectors, it does have Killing vectors 
which are timelike in half of the spacetime and exhibit a bifurcate Killing horizon structure. One 
can introduce static coordinates which cover each one of the two patches where the Killing vec- 
tor is timelike (hence one quarter of the full de Sitter space) and correspond to the following 
parameterization of the hyperboloid: 



X° = sinh(t), 
X' = -roj\ i = 1,2,3, 

= -Vl -r2 cosh(t), (66) 

where < r < 1 and w*, which satisfy the condition 5ijUJ^uj^ = 1 and define a 2-sphere, can be 
parametrized as 

oj^ = sin(6'2) sin(0i), 
= sin(6'2) cos(0i), 
^3 = cos(02). (67) 

The line element in these coordinates is given by 

ds^ = -{l-r^)dt^ + j:^^^ + r^dnl (68) 

where dil2 = ^^2 + sm'^{62) d0f. Four independent isometry generators are manifest in these 
coordinates: three associated with the 0(3) group of the 2-sphere and one corresponding to the 
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time-translation symmetry t — )• t + a. The worldlines of constant r are uniformly accelerated 
timelike curves and only r = is a geodesic (with proper time t), which coincides with the 
same geodesic considered in the previous two subsections; r = 1 corresponds to the event horizon 
surrounding this geodesic observer. 

Note that there is an analogous static patch, antipodal to the one considered so far, which 
can be covered by the same kind of coordinates by changing the sign of both and X'^ in the 
above parameterization. Furthermore, the remaining two quadrants of de Sitter can be covered by 
extending the range of the radial coordinate to 1 < r < oo and making the following changes in 
the embedding parameterization described by eq. ([66]) : 

X° = ±\/r2 - 1 cosh(t), 

X'^ = =FVr2 _ I sinh(t). (69) 

The two additional quadrants covered by this parameterizations are not static since the timelike 
and spacelike character of dt and dr are interchanged. 

We conclude this appendix briefly describing the relation between the static coordinates (as 
well as their extensions) and the coordinate systems of the previous subsections. 

a. Relation to spatially flat coordinates 

Comparing the expressions for X* in eqs. ()63p and ()66p one immediately finds 

I 

r = = a{r]) \x\. (70) 

V 

Similarly, by comparing the expressions for X^ — X^, one gets 

- i = e* ^/l-r^ or e"^* = - I^P- (71) 
V 

Taking the logarithm, the following relation between the comoving time of spatially flat coordinates 
and the static time is obtained:^ 

T = t + ^ln(l-r2). (72) 

They both coincide along the worldline at r = and correspond to the proper time of that geodesic. 
Furthermore, from these relations it is straightforward to see that the isometry transformation 
defined in eqs. (|65p corresponds to the time translation symmetry for t ^ t + a keeping r fixed in 
static coordinates. 



■* If one extends the {r, t} coordinates outside the static patch using eqs. H69|l , there the relation becomes r = 
-t+(l/2) ln(r2-l). 
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b. Relation to global coordinates 



Finally, comparing (with i = 1,2,3) as well as / X'^ in eqs. (j6ip and (j66p . one finds the 
following relation between global and static coordinates: 

r = cosh(r)sin(03), (73) 

- («) 

If one extends the static coordinates outside the static patch, the last relation becomes 

, , , tanh(T) 

coth t = 75 

cos(6'3) 



APPENDIX B: REGULAR EXTENSION OF THE STATIC SOLUTION 

The static solution for a free-falling charge in de Sitter found in Sec. |V] is directly related to a 
solution of the Euclidean version of the Klein-Gordon equation on the 4-sphere. This fact can be 
exploited to show that for > the solution regular on the horizon can always be matched to 
the solution with a source at r = 0. This will be shown in the first subsection of this appendix by 
arguing that if the contrary were true, there would be a solution regular everywhere on the sphere, 
and then proving that this is impossible for > 0. In the second subsection we will discuss the 
extension of this solution to the full de Sitter spacetime and its regularity. 



1. Equivalent problem on the 4-sphere 



A 4-sphere can be obtained by analytic continuation from de Sitter space by taking t = i9 in 
eq. (j68|) and requiring 6 to have a period of 27r so that there is no conical singularity at r = 1. This 
becomes perhaps even clearer if we introduce the change of coordinates r = sin(x), with < x ^ "^j 
so that the line element becomes 

ds^ = cos"^ {x)d9^ + dx^ + sm'^{x)dnl (76) 

One can proceed in the same way to obtain the Euclidean version of the Klein-Gordon equation (|36p . 
In the new variables the source term is given by {q/'^'^X^)^{x) + {q/^'^{x ~ ~ ^)) where 

the second term would correspond to the second static patch mentioned in Appendix A. 3 and is 
required here for continuity in 6. If one looks for solutions of the form (f){x) = 'A(^) (iii fact, for 
> they are necessarily of this form), the equation satisfied becomes equivalent to eq. (j37p .^ 
Thus, one can make use of the results in eqs. (j38p -()4ip for the expansions around the regular 
singular points at x = and x = 7r/2 (corresponding respectively to r = and r = 1). Taking into 



^ Note that it makes sense to identify 4>{x) = only if 4>{x) = 4>{'^ ~ x)- This is possible because the smooth 
function that we will find below is an even function of (x — 7r/2). Hence, from now on we will focus on x < ""/S 
and use (pix) = ~ x) for X ^ 7r/2. 
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account that r = sin(x), it is clear that (j)i{x) from eq. ()38p is an even function of x with no spike 
at X = 0, where it is smooth. Similarly, taking into account that 1 — r = 1 — cos (x — '?i"/2), one 
can see that 4>i{x) from eq. (I40p is an even function of (x — 7r/2) with no spike at x = ^/"^^ where 
it is smooth. 

The advantage of working on the 4-sphere is that we can easily prove that for M2 > 0, CMx) 
always matches a linear combination A(j)i{x)~^B(j)2{x) with non-vanishing B. Indeed, if B vanished, 
C(j)i{x) would match and would correspond to a homogeneous solution of eq. (i37|l which 

would be regular everywhere on the 4-sphere, but one can show that this is impossible as follows. 
The homogeneous solutions of eq. (|37p are solutions of the Euclidean version of the Klein-Gordon 
equation: 

(A4 - = 0, (77) 

where A4 is the Laplace operator on S^, whose eigenvalues are —L{L + 3) with L = 0, 1, . . . 
Therefore, one cannot have a regular solution for M"^ > since that would correspond to an 
eigenfunction of the Laplacian with a positive eigenvalue. In contrast, for M"^ = one does have 
the obvious solution (p = const., which explains why our argument does not work for the massless 
minimally coupled case and why one cannot find a static solution regular on the horizon associated 
with a freely-falling charge in that case. 



2. Extension to full de Sitter spacetime 



One can obtain de Sitter space by analytic continuation from S"^ reversing the procedure at 
the beginning of the previous subsection. In that way, one is naturally led to a static patch of de 
Sitter, but one can go to full de Sitter by extending the static coordinates to r > 1, as explained in 
Appendix A. 3. Alternatively, one can consider an analytic continuation of the form iT = 64 — Tr/2 
to the global coordinates. On the 4-sphere the source is located on an equatorial worldline with 
63 = {0, vr} and arbitrary values of 64, which becomes two disconnected antipodal worldlines in de 
Sitter located respectively at 6*3 = and 6*3 = vr, and arbitrary values of T. 

Since we found a function which is regular everywhere on (outside the source), the analytic 
continuation to full de Sitter space should be regular too. This can be explicitly checked by using 
eq. (|73p to express in terms of global coordinates the solution (j)s{r) found in Sec. IV Ai The only 
potentially problematic value is r = 1 and, as seen in eq. (j40p . the solution around that value can 
be expanded as a convergent series in powers of (1 — r), which can be written in global coordinates 
as 

1 1 \./rr\ ■ fa \ 1 sin((93) cos(6l3 - 7r/2) 

1 - r = 1 - cosh(r) sm(6'3) = 1 — - = 1 — , (78) 

cos(A) cos(A) 

where the neighborhood of r = 1 can be more easily studied in terms of the conformal time A, 
since r = 1 then corresponds to ^3 — tt/2 = ibA mod vr. The solution is, thus, manifestly smooth 
everywhere outside the sources. Note that in order not to have a spike at ^3 — 7r/2 = A = it is 
important that (psi^) is an even function of ^3 — 7r/2 and A, which follows from eq. (|73p and the 
fact that the solution in eq. pUj) is analytic in (1 — r). 
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Besides the neighborhood of r = 1 one can also study the behavior at late times (and fixed 
comoving distance) or large physical distance from the sources (and fixed time) by considering the 
solution (psif) for large values of r. More specifically, one can introduce the change of variables 
w = 1/r in eq. ()37p (note that for large r the source is irrelevant) and see that u) = is a regular 
singular point. One can, therefore, make use of the Frobenius method and Fuchs's theorem to 
express the solutions as a series expansion in powers of 1/r, analogously to what was done in 
Sec. IV A[ The two possible solutions are given by 

Mr) = r-^- (! + •••)' (79) 
Mr) = r-^+ (! + •••)' (80) 

where X± = 3/2 it y^9/4 — (notice that for > their real part is always positive) and 
the remaining terms involve positive integer powers of 1/r times a possible Inr factor in some 
cases. Using r = cosh(T) sin(03) one can straightforwardly obtain the late-time and large-distance 
behaviors in global coordinates. (The late-time behavior coincides with that obtained in Sec. IV Bl 
for the homogeneous solutions of the Klein-Gordon equation because the source played no role here.) 
The extension of the static solution will correspond to a unique linear combination D 4>i{r)+E 4>2{r) 
whose coefficients should be determined by matching it to the regular solution near the horizon, 
given by eq. (|40|) . with the right normalization factor. This would require solving exactly (|37p . 
Lacking an exact solution, we can still expect that the large r behavior will be determined by 
(/)i(r), which dominates over i?i>2(r) at large r, as long as L> / (that this is the case could perhaps 
be shown using approximate methods to match the solutions). For the particular example of a 
conformal field (corresponding to = 2), we know that the exact solution is proportional to 1/r 
and D / indeed (while E = {)). 
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